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Abstract 

This  paper  develops  a  model  for  the  study  of  solid  oxide  fuel  cell  (SOFC)  electrical  efficiency.  The  chemically  reacting  flow  model 
considers  the  coupled  behavior  of  fluid  mechanics,  thermodynamic  equilibrium,  and  transport  of  electroactive  species  through  the  elec¬ 
trolyte.  Fuel  cells  are  treated  as  plug  flow  tubular  reactors  with  no  upstream  diffusion.  Minimization  of  free  energy  is  used  to  calculate 
gas  composition  along  the  length  of  the  fuel  channel  without  specifying  a  set  of  reactions.  The  model  serves  to  place  an  upper  bound  on 
electrical  efficiencies.  The  usefulness  of  the  model  is  demonstrated  by  simulating  SOFC  electrical  efficiency  as  a  function  of  operating 
voltage  for  in  situ  reforming  of  four  different  fuel  mixtures,  including  steam  reforming  and  partial  oxidation  of  CH4,  syngas,  and  pure 
CH4. 
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1.  Introduction 

The  commercial  feasibility  of  fuel  cells  rests  on  the  cost 
of  the  fuel  cell  system  and  operating  cost.  The  latter  is  deter¬ 
mined  mostly  by  electrical  efficiency  and  fuel  cost.  Anode- 
supported  solid  oxide  fuel  cells  (SOFCs)  are  especially  attrac¬ 
tive  since  they  are  both  fuel-flexible  and  potentially  highly 
efficient  [1]. 

Our  paper  reports  the  limiting  electrical  efficiency  of 
SOFCs  as  constrained  by  thermodynamics  for  three  dif¬ 
ferent  fuel  mixtures  and  pure  CFLp  In  situ  processing  of 
syngas  (66.7%  FL  and  33.3%  CO),  CH4  by  steam  re¬ 
forming,  and  CH4  by  partial  oxidation  are  reported.  POX 
is  especially  attractive  since  it  obviates  the  need  to  inject 
steam  into  the  fuel  stream.  The  model  developed  here  de¬ 
termines  the  maximum  intrinsic  electrical  efficiency  ob¬ 
tainable  in  a  SOFC  for  any  hydrocarbon  fuel  mixture. 
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The  essential  question  addressed  is  whether  or  not  ce¬ 
ramic  fuel  cells  can  reach  electrical  efficiencies  higher  than 
about  60%,  in  order  to  be  competitive  with  utility  power 

[2,3]. 

Performance  of  the  fuel  cell  is  modeled  on  the  basis  of  gas 
composition  in  the  anode  channel.  Although  the  gas  com¬ 
position  at  the  cathode  is  also  important,  flow  in  the  cath¬ 
ode  channel  (generally  ambient  air)  is  typically  many  times 
greater  than  that  in  the  anode  channel.  Thus,  gas  composi¬ 
tion  along  the  cathode  channel  may  be  treated  as  constant. 
With  fast  chemical  kinetics  and  minimal  mass  transport  lim¬ 
itations,  gas  in  the  anode  channel  will  approach  the  com¬ 
position  predicted  by  thermodynamic  equilibrium.  There¬ 
fore,  ionic  transport  of  the  ions  (O2-)  through  the  electrodes 
and  electrolyte  becomes  the  rate-determining  step  in  our 
model. 

An  advantage  of  the  minimization  of  free  energy  method 
[4]  is  that  it  only  requires  the  free  energy  of  reactant  and 
product  species.  Knowledge  of  reaction  mechanisms  or  equi¬ 
librium  constants  is  not  required,  a  priori.  Nonetheless,  all 
major  species  must  be  accounted  for.  Sasaki  and  Teraoka 
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[5]  examined  a  wide  range  of  possible  SOFC  fuels  and 
found  five  primary  fuel  cell  gases,  H2,  H2O,  CO,  CO2, 
CH4,  as  well  as  graphite.  Moreover,  they  predicted  the  open 
circuit  potential  of  the  thermodynamic  equilibrium  mix¬ 
ture  in  terms  of  the  C-H-0  ratio  and  used  Gibbs  min¬ 
imization  to  calculate  equilibrium  gas  mixtures.  In  a  ki¬ 
netic  study  of  steam  reforming  of  methane  and  water  gas 
shift  (WGS)  reactions  on  a  Ni/MgAl204  catalyst,  Xu  et 
al.  [6]  proposed  an  11 -reaction  mechanism,  which  included 
the  same  five  gas-phase  species  reported  by  Sasaki  and 
Teraoka. 

Demin  et  al.  [7,8]  developed  a  thermodynamic  model 
for  estimating  the  system  efficiency,  including  electrical  ef¬ 
ficiency,  of  SOFCs  operating  on  steam-methane  mixtures. 
Their  model  developed  an  algorithm  for  computing  the  gas 
composition  along  the  anode  channel.  The  algorithm  assumes 
the  same  gas  phase  species  as  Sasaki  and  Teraoka  and  Xu  et 
al.  and,  unlike  our  model,  assumes  two  reactions:  steam  re¬ 
forming  reaction  and  water  gas  shift.  An  important  distinction 
between  the  Demin  et  al.  models  and  our  model  is  that  they 
use  the  average  EMF  to  compute  electrical  efficiency  and  as¬ 
sume  that  all  processes  occur  under  reversible  conditions,  i.e. 
no  load.  In  contrast,  our  model  calculates  the  power  gener¬ 
ated  at  each  point  in  the  cell  from  the  operating  voltage  and 
local  current.  We  assume  the  zero  current  condition  only  at 
the  end  of  the  cell.  Since  our  model  includes  ohmic  losses, 
electrical  efficiencies  reported  in  this  paper  are,  therefore, 
more  realistic. 

Accounting  for  high  mass  transport  resistance  in  the  an¬ 
ode  support  and  slow  surface  chemistry  requires  a  more  com¬ 
plex  model.  Moreover,  inclusion  of  these  effects  also  requires 
detailed  knowledge  of  the  porosity  and  permeability  of  the 
anode  support  together  with  well-developed  surface  reaction 
mechanisms.  Zhu  et  al.  developed  a  SOFC  model  which  in¬ 
cludes  both  channel  flow  and  porous  media  flow  coupled  with 
heterogeneous  surface  chemistry  in  the  anode  support  [9,10]. 
In  Fig.  9  of  the  paper,  Zhu  et  al.  [10]  compared  gas  composi¬ 
tion  and  current  densities  along  the  anode  channel  assuming 
either  heterogeneous  kinetics  or  chemical  equilibrium.  In  the 
range  of  interest  where  commercial  fuel  cells  are  operated, 
that  is,  at  high  thermodynamic  and  Faradaic  efficiency,  both 
kinetic  and  equilibrium  models  predict  nearly  equivalent  fuel 
cell  performance. 

A  difficulty  with  kinetic  models  is  that  the  validity  of  the 
predictions  depends  on  detailed  a  priori  knowledge  of  rate 
constants.  Such  knowledge  is  often  not  available  for  the  con¬ 
ditions  of  interest.  In  the  present  model,  gas  in  the  anode 
channel  is  assumed  to  be  at  chemical  equilibrium  along  the 
entire  length  of  the  fuel  cell.  Moreover,  by  providing  an  ar¬ 
bitrarily  long  fuel  cell,  the  fuel  cell  is  able  to  reach  reversible 
conditions  (zero  current)  by  the  time  the  fuel  exits  the  cell. 
We  do  not  seek  to  replace  kinetic  models.  Instead,  our  model 
serves  to  establish  an  upper  bound  of  electrical  efficiency 
based  on  the  Second  Law  rather  than  complex  overpoten¬ 
tials  which  are  fuel  cell  specific  and,  in  many  cases,  poorly 
understood. 


2.  Model  development 

Composition  of  the  gas  flow  in  the  anode  channel  de¬ 
pends  on  the  coupled  behavior  of  fluid  mechanics,  thermo¬ 
dynamic  equilibrium,  and  transport  of  electroactive  species 
through  the  electrolyte.  Gases  in  the  anode  channel  are 
treated  as  a  plug  flow  (radially  well-mixed  with  no  ax¬ 
ial  diffusion).  For  the  very  low  flow  rates  considered  here 
(velocities  of  the  order  of  3cms-1  and  Reynolds  num¬ 
ber  10),  pressure  losses  are  insignificant  over  the  length 
of  the  tube.  Gas  pressure  and  temperature  in  both  anode 
and  cathode  channels  is  assumed  to  be  isobaric  and  isother¬ 
mal. 

Fig.  1  illustrates  two  arrangements  for  anode- supported 
fuel  cells.  Fuel  mixtures  flow  through  the  interior  anode  chan¬ 
nels.  Air  flows  around  the  outside  of  the  tube  (left)  or  through 
the  bottom  channel  (right).  On  the  left  side  of  Fig.  1,  the  chan¬ 
nel  is  cylindrical  and  the  electrolyte  extends  around  the  entire 
perimeter.  For  the  cylindrical  tubes,  the  fuel  cell  stack  con¬ 
sists  of  multiple  tubes  connected  to  a  fuel  manifold.  For  the 
planar  configuration  (right  side  Fig.  1),  three  of  the  walls  are 
metal  interconnects. 

The  anode  support  and  electrolyte  boundary  occupy  the 
bottom  of  the  anode  channel.  Although  our  model  is  devel¬ 
oped  for  tubular  geometry,  it  is  general,  and  may  be  applied 
to  other  configurations  where  the  cross-section  perpendicular 
to  flow  and  electrolyte  area  are  known. 

The  anode  support,  typically  Ni/8YSZ  cermet,  is  a  good 
electrical  conductor  and  each  tube  operates  at  a  single  voltage, 
Well  • 

Transport  of  O2-  ions  from  cathode  to  anode  alters  the 
local  gas  composition  along  the  anode  channel  by  oxidiz¬ 
ing  the  fuel.  For  a  sufficiently  long  tube,  the  extent  of  fuel 
oxidation  depends  on  the  cell  operating  voltage,  Vcen.  As 
the  mixture  becomes  oxidized,  the  reversible  voltage,  Ev , 
drops  until  it  reaches  VCq\\.  Beyond  this  point,  additional 
ionic  current  cannot  flow,  and  it  is  not  possible  to  oxidize 
the  fuel  any  further.  Faraday’s  Law  relates  the  local  molar 


Tubular  Planar 


Fig.  1.  Anode- supported  SOFCs.  Left,  a  tubular  arrangement  where  the 
entire  tubular  anode  support  is  covered  with  electrolyte;  right,  a  square 
channel  arrangement  where  one  side  of  the  channel  is  covered  with 
electrolyte. 
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Transport  of  species  in  the  anode  channel 

j-1  j-1/2  j  j-h  1/2  j+1 

— • - 1 - • - 1 - • - 

x  x  +  Ax 


Fig.  2.  The  fuel  cell  is  modeled  as  a  plug  flow  reactor  with  no  upstream 
diffusion.  The  gas  entering  each  volume  element  reaches  chemical  equilib¬ 
rium. 


flow  of  O2  ions  through  the  electrolyte  to  the  local  electric 
current, 


i  =  z*F—,  (1) 

at 

where  i  is  in  A,  ze  =  —2  the  charge  transfer  per  ion,  F  Fara¬ 
days  constant,  and  n  is  moles  of  O2-  ions. 

Our  model  uses  a  shell  balance  to  segment  the  tube  into 
discrete  volume  elements.  Fig.  2  represents  the  shell-balance 
control  volumes  of  elements  j  and  j  +  1 .  Three  control  sur¬ 
faces  account  for  transport  of  species  into  each  control  vol¬ 
ume.  A  mole  balance  of  species  i  entering,  leaving,  reacting, 
and/or  accumulating  within  a  control  volume  is  shown  in  Fig. 
2,  where  hij-1/2  is  the  molar  flow  of  species  i  into  the  el¬ 
ement  7,  Ri  j  the  rate  of  creation  or  annihilation  of  species 
i  caused  by  chemical  reactions  in  element  j,  ho2  .  the  molar 

flow  of  O2  due  to  transport  of  O2-  ions  through  the  elec¬ 
trolyte,  and  hij+1/2  is  the  molar  flow  out  of  element  7.  For 
the  steady  state  case,  the  net  rate  of  accumulation,  driij/dt ,  is 
zero  for  all  species. 

Summing  over  all  the  species,  A,  in  element  7,  at  steady 
state,  the  molar  flows  in  Fig.  2  become, 

N  N  N 

ni, j-1/2  +  no2j  —  y^/ij,7+i/2  +  ^  Ri,j  —  0.  (2) 

i=l  i= 1  i=\ 


ture  flow, 

N  N  N 

r,  Rij — uj — y^^,7-y2>  (3) 

i=\  i=  1  i=  1 

where  positive  Ri  j  represents  creation  of  species  i  and  neg¬ 
ative  Rf  j  represents  annihilation  of  species  i.  The  rate  of 
fuel  oxidation  depends  on  the  rate  of  transport  of  O2-  ions 
through  the  electrolyte.  Nevertheless,  we  can  immediately 
find  the  equilibrium  composition  in  each  element  j  as  a  func¬ 
tion  of  the  amount  of  O2  added  to  the  anode  channel.  Our 
model  relies  on  HSC  software  [12]  to  minimize  the  Gibbs 
free  energy  in  each  volume  element  (Eq.  (3)).  Combining 
Eqs.  (2)  and  (3)  yields, 

N  N 

y^^eq  ij  +  no2J  =  y>7+l/2,  (4) 

i=l  i= 1 

which  is  the  non-equilibrium  molar  flow  leaving  volume  ele¬ 
ment  j  and  entering  element  j  +  1.  Thus,  the  calculation  pro¬ 
ceeds  in  a  marching  fashion  from  the  channel  inlet,  j  =  1 ,  to 
the  last  volume  element  of  the  cell,  j  =  L. 

Generally,  moles  are  not  conserved  in  the  control  volume 
due  to  chemical  reactions  and  net  addition  of  O2  through 
the  electrolyte.  As  a  result,  gas  velocity  increases  along  the 
length  of  the  tube.  For  an  ideal  gas,  the  velocity  of  the  gas 
into  element  j  is  stated  as, 

(Pj_i/2)(cms-1)  =  22414(cm3  mol-1) 

T  P()  E^i'Vj-i/ataols-') 

X  273  P  A  (cm2) 

(5) 


where  T  is  temperature  (K),  P  pressure,  fiij-1/2  molar  flow 
of  species  i  through  the  control  surface  at  7  —  1  /2,  A  the  total 
number  of  species,  and  A  is  the  cross-sectional  area  of  the 
channel. 

In  each  element 7,  the  equilibrium  gas  composition  is  per¬ 
turbed  by  the  transport  of  O2-  ions  through  the  electrolyte. 
Combining  Ohms  Law  and  Eq.  (1),  the  current  produced  by 
a  SOFC  volume  element  is, 


Eq.  (2)  is  a  form  of  the  general  mole  balance  equation, 
which  is  described  elsewhere  [11].  We  assume  that  the  molar 
flow  into  each  volume  element,  hij- 1/2,  immediately  reaches 
equilibrium  upon  entering  the  control  volume.  Chemical  re¬ 
actions  create  and/or  destroy  species  to  reach  equilibrium 
composition,  which  is  calculated  by  Gibbs  minimization. 
Chemical  equilibrium  itself  depends  upon  gas  composition, 
temperature  T,  and  pressure  P.  To  reach  chemical  equilib¬ 
rium,  the  rate  at  which  each  specie  is  created  or  destroyed 
must  be  equal  to  the  difference  between  the  non-equilibrium 
flow  coming  into  the  control  volume  and  the  equilibrium  mix¬ 


where  Vceii  is  the  cell  voltage,  2  x  ze  =  4  the  charge  transfer 
for  a  O2  molecule,  R0hmj  the  apparent  electrical  resistance 
of  the  element,  F  Faradays  constant,  and  fio2j  is  the  molar 
flow  of  O2.  The  reversible  voltage  for  element  j  in  Eq.  (6)  is, 


Erj  =  E°j  + 


RT 
4 ~F 


ln(Fo?(c))  +  —  In 


RT  (  Pn2j(a) 


2  F  \Pn2o,j(a ) 


where  c  and  a  indicate  partial  pressures  at  the  anode  and  cath¬ 
ode,  respectively.  The  reversible  voltage  of  the  cell  depends 
on  the  equilibrium  potential  between  the  electrodes,  and  is 
taken  as  the  open  circuit  voltage  of  element  7.  The  operat- 
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in g  voltage  of  the  cell,  Vceu,  must  be  less  than  the  reversible 
voltage,  7+,  for  current  to  flow  spontaneously. 

As  an  electrical  load  is  applied  to  a  fuel  cell,  the  voltage 
across  the  electrolyte  can  be  represented  most  generally  as, 


^cell  —  Ev  7?conc,a  Osl  Oohm  Oconc,c  Oc 

“^interface  —  ^leakage?  (8) 

where  77COnc,a  and  rjc onc,c  are  the  concentration  overpotentials 
at  the  anode  and  cathode  due  to  diffusion  resistance,  rja  and 
Oc  the  activation  polarizations  at  the  anode  and  cathode,  ijohm 
the  ohmic  polarization,  interface  represents  surface  polariza¬ 
tion  phenomena,  and  and  leakage  is  the  polarization  loss  due 
to  electron  leakage  [9].  For  the  model  presented  here,  we 
lump  all  of  the  linear  polarization  phenomena  into  the  ohmic 
polarization  term  r/ohm  and  discard  all  non-linear  polarization 
effects.  Activation  polarizations,  which  do  not  depend  on  cur¬ 
rent,  are  assumed  to  be  much  smaller  than  ijohm,  and  are  not 
considered.  Consequently,  the  current-voltage  (IV)  charac¬ 
teristic  curve  is  treated  as  a  straight  line,  and  each  R0hm,  j> 
in  units  of  ohms,  is  constant  within  each  element  j.  For  the 
isothermal  conditions  investigated  here,  the  area-specific  re¬ 
sistance  (ASR),  in  units  of  Q  cm2,  is  constant  for  all  volume 
elements. 

Fig.  3  illustrates  IV  curves  for  two  adjacent  volume  el¬ 
ements  j  and  y'+l.  Each  element  j  has  its  own  reversible 
voltage  determined  from  the  gas  composition  in  the  volume 
element.  The  upper  curves  represent  element  j,  which  has  a 
higher  reversible  voltage  Ev.  A  short  distance  down  the  chan¬ 
nel,  element  j  +  1 ,  the  reversible  voltage  is  lower  and  the  IV 
curve  is  shifted. 

The  excess  potential  A  V)  is  equal  to  the  difference  be¬ 
tween  the  reversible  voltage,  Etj,  and  the  cell  voltage,  Vceu. 
Ohms  Law  holds  for  each  volume  element  and  the  excess 
voltage,  A  V,  is  related  to  the  current  ij  by, 


where  ij  and  R0 hm,;  (units  of  ohms)  are  the  current  and  elec¬ 
trical  resistance  of  volume  element  7,  respectively.  The  re¬ 
sistance  in  each  element,  which  depends  on  the  ASR  and 
electrolyte  area  contacting  the  gas  in  the  anode  channel,  is 
stated  as, 


ASR 

TtDAxj  ’ 


(10) 


where  ttDAxj  is  the  wetted  electrolyte  area  A  j  (see  Fig.  1). 

Both  top  and  bottom  panels  in  Fig.  3  show  ohmic  losses 
in  each  element  as  a  dark  shaded  triangular  region  above 
the  cell  voltage  Vceii.  The  power  in  each  element  j  is  shown 
as  square  shaded  regions  below  the  operating  voltage  Vceii. 
Importantly,  Fig.  3  shows  that  the  electrical  power  of  every 
element  j  includes  all  ohmic  losses. 

Consider  two  alternatives  to  calculate  the  molar  flow 
of  O2,  no2J>  into  ea°h  volume  element.  We  could  fix 
Axj,  (making  the  volume  and  electrical  resistance  of  each 
segment  the  same,  see  Fig.  2)  and  solve  for  the  current 


Fig.  3.  Each  segment  along  the  length  of  the  fuel  cell  tube  has  its  own  TV 
curve;  curves  for  j  and  j  +  1  are  shown  here.  In  the  top  panel,  the  length  of 
each  volume  element  Ax  j  and  resistance  is  held  constant  for  all  elements. 
The  lower  panel  shows  the  IV  curves  for  adjacent  elements  with  constant 
current.  In  both  top  and  bottom  panels,  the  ASR  is  constant. 


by  means  of  Eq.  (6).  Alternatively,  we  could  add  a  con¬ 
stant  amount  of  O2  to  each  volume  element,  and  allow 
the  segment  length  to  vary.  The  next  sections  discuss  each 
approach. 

2.7.  Dimensional  solution 

The  current  and  molar  flow  of  O2-  ions  in  each  volume  el¬ 
ement  will  vary  with  the  excess  potential,  A  V),  while  the  elec¬ 
trical  resistance  remains  fixed.  The  solution  is  then  achieved 
in  a  marching  fashion  from  the  inlet  to  the  outlet  in  steps  of 
Ax.  Importantly,  this  approach  requires  fuel  cell  specific  ge¬ 
ometrical  information,  tube  diameter,  ASR,  electrolyte  area 
(through  Eq.  (10)),  and  operating  conditions. 

As  shown  in  the  top  panel  of  Fig.  3,  the  slopes  of  the  j 
and  7  +  I  curves  are  equal  since  the  electrical  resistance  in 
each  element  is  constant.  As  fuel  in  the  anode  channel  is  con¬ 
sumed,  a  family  of  parallel  7V  curves  is  formed,  each  shifting 
downwards  as  the  reversible  voltage,  Ev ,  drops.  Current  in 
each  segment  is  found  from  the  intersection  of  Vceii  and  the 
IV curve.  Ideally,  Erj  =  Vceii  in  the  final  element  (7  =  L)  and 
the  current  is  zero.  In  practice,  the  reversible  voltage  drops 
very  rapidly  as  the  fuel  nears  complete  oxidation  and  the  final 
element  is  identified  as  last  element  having  positive  excess 
voltage  AVj. 
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2.2.  Dimensionless  solution 


The  total  channel  length  is  found  by, 


In  the  dimensionless  approach,  the  molar  flow  of  O2  into 
each  element  is  fixed  and  we  solve  for  the  reversible  voltage 
Erj  of  each  element  j.  From  Faradays  Law,  Eq.  (1),  it  can 
be  seen  that  current,  ij,  and  molar  flow  of  oxygen,  ho2j,  are 
proportional. 

In  the  bottom  panel  in  Fig.  3,  the  molar  flow  of  O2  into 
each  element,  no2,j,  is  fixed.  The  length  of  each  element, 
Ax  j,  changes  to  adjust  the  element  resistance  and  maintain 
a  fixed  current.  As  fuel  is  consumed  the  reversible  potential, 
Fr,  drops  and  the  IV  curve  rotates.  A  family  of  IV  curves, 
which  pivot  about  a  single  point,  are  formed.  Ideally,  in  the 
last  element,  the  reversible  voltage  Fr  reaches  the  cell  voltage 
Vceii  and  the  resistance  reaches  zero.  Numerically,  the  end  of 
the  cell  is  reached  when  a  discrete  addition  of  O2  results  in 
a  negative  excess  potential  A  V  and  the  slope  of  the  IV  curve 
becomes  positive. 

The  equilibrium  gas  composition  and  reversible  voltage 
along  the  fuel  cell  is  known  as  a  function  of  O2  added  per 
mole  of  inlet  fuel  mixture,  i.e.  the  extent  of  the  oxidation  re¬ 
action.  Moreover,  the  stoichiometric  amount  of  O2,  needed  to 
fully  oxidize  the  inlet  mixture,  is  also  known.  Therefore,  the 
problem  can  be  scaled  on  a  molar  basis,  with  the  amount  of 
O2  added  to  the  inlet  mixture  ranging  from  zero  to  stoichio¬ 
metric  (as  a  practical  matter,  accuracy  of  the  solution  may 
be  improved  by  adapting  the  step  size  of  ho2j  to  the  rate  of 
change  of  the  reversible  voltage  Fr).  For  a  given  operating 
voltage,  Vceii,  both  the  excess  potential,  A  V ,  and  total  current 
are  also  known.  With  reversible  voltage,  total  current  and  fuel 
flow  in  hand,  the  electrical  efficiency  is  easily  calculated,  and 
is  defined  as, 


«p  = 


Aotal 


Ef=  1  ntAHi 


(11) 


where  Aotai  is  the  total  electrical  power,  hi  the  inlet  molar 
flow  rate  of  fuel  species  i  per  mole  of  total  gas  flow  at  the 
inlet,  A  the  number  of  fuel  species,  and  A  Hi  is  the  enthalpy 
for  complete  oxidation  of  species  i.  The  operating  voltage  is 
constant,  the  electrical  power  from  a  single  volume  element 

j  is, 


Pj  =  L' Vceii-  (12) 

Therefore,  the  total  electrical  power  for  the  fuel  cell  is, 

L 

Aotal  =  4 F Vcei)  J>o2,„  (12) 

7=1 


where  F  is  Faradays  constant,  ho2j  the  molar  flow  of  O2 
into  each  element  j ,  and  L  is  the  last  volume  element  with  a 
positive  excess  potential  AV.  The  electrical  efficiency  (Eq. 
(11))  becomes, 


4 F Vceii  X] 7=1  ^02,7 
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where  L  is  the  number  of  volume  elements  and  Ax  j  is  the 
length  of  element  j.  Of  course,  Axj  goes  to  infinity  as  the  ex¬ 
cess  potential  AV  j  and  current  density  approach  zero.  There¬ 
fore,  a  small  A  V  must  be  chosen  to  truncate  the  solution. 

Since  the  electrical  efficiency  of  the  fuel  cell  depends  only 
on  temperature,  composition,  and  operating  voltage,  a  more 
general,  dimensionless,  solution  to  the  problem  is  obtained. 
The  general  solution  may  then  be  used  to  calculate  the  length 
of  each  volume  element  based  on  fuel  cell  specific  informa¬ 
tion,  i.e.  ASR,  channel  cross-section,  electrolyte  area  per  unit 
length,  and  molar  flow  at  the  inlet. 

The  current  per  mole  of  inlet  flow  of  each  element  j  is  used 
to  scale  the  solution  to  the  specific  geometry  (see  Fig.  1)  of  a 
particular  fuel  cell.  For  a  fixed  value  ho2j,  the  length  of  the 
volume  element,  Ax,  changes  to  accommodate  the  varying 
current  density.  Thus,  the  length  of  each  volume  element  Ax  j 
is  needed  to  scale  the  solution.  Both  the  molar  flow  of  O2  and 
excess  voltage  are  known  and  the  length  of  volume  element 
j  is  found  as, 


ASR  4Ffio2,7 

n D  ( Erj  —  Veen) 


(16) 


Eq.  (16)  together  with  Erj  is  used  to  calculate  current  den¬ 
sities  and  power  profiles.  From  the  equilibrium  gas  compo¬ 
sition  of  each  element  j ,  one  may  obtain  mole  fraction  and 
velocity  profiles  inside  the  anode  channel  (Eqs.  (6)  and  (16)). 


3.  Results  and  discussion 

Our  model  can  be  applied  to  a  wide  range  of  fuel  cell 
conditions  to  provide  useful  insight  into  various  design  fuel¬ 
ing  options.  In  the  following  example,  different  methods  of 
processing  methane  are  considered. 

Fig.  4  shows  the  electrical  efficiency  as  a  function  of  oper¬ 
ating  voltage  for  four  different  inlet  mixtures  (composition  in 
mole  fractions),  pure  CH4,  POX  (0.408  CH4, 0. 123  O2, 0.059 
H2O,  and  balance  N2),  syngas  (0.667  H2  and  0.333  CO),  and 
steam  reforming  (0.667  H2O  and  0.333  CH4).  Syngas  com¬ 
position  varies  depending  on  the  hydrocarbon  source.  The  2: 1 
H2:CO  ratio  represents  the  inlet  fuel  composition  one  would 
expect  from  steam  reforming  a  50/50  mixture  CH4  and  coal 
at  high  temperature.  In  all  cases,  temperature  is  isothermal 
and  maintained  at  800  °C. 

Pure  CH4  (triangles)  is  shown  to  give  a  bounding  case  and 
results  in  the  highest  electrical  efficiency,  75%.  In  practice, 
pure  CH4  cannot  be  used  because  of  coke  formation.  The 
maximum  electrical  efficiency  for  the  partial  oxidation  case, 
63%  at  8 14  mV,  is  considerably  lower  than  that  for  steam 
reforming,  72.7%  at  797  mV  (shown  with  squares).  Syngas 
(shown  with  crosses)  yields  the  lowest  electrical  efficiency, 
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Cell  operating  voltage  (V) 

Fig.  4.  SOFC  electrical  efficiency  as  a  function  operating  voltage  for  four 
different  gas  mixtures.  Temperature  is  maintained  at  800  °C  and  pressure  is 
1  bar. 

57%  at  815  mV.  As  shown,  steam  reforming  results  in  a  lower 
open  circuit  voltage  and  higher  electrical  efficiency  than  ei¬ 
ther  syngas  or  POX.  An  important  finding  from  the  model  is 
that  both  internal  steam  reforming  and  POX  can  be  employed 
in  fuel  cells  with  >60%  electrical  efficiency  while  the  use  of 
2:1  syngas  cannot. 

Consider  the  low  efficiency  regions  in  Fig.  4.  For  op¬ 
erating  voltages  ranging  from  zero  to  around  700  mV  effi¬ 
ciency  responds  linearly  with  the  operating  voltage  because 
oxidation  of  the  fuel  is  nearly  complete.  This  is  a  region 
characterized  by  high  fuel  efficiency  (Faradaic  efficiency) 
but  low  thermodynamic  efficiency.  A  second  region  of  low 
efficiency  occurs  at  high  operating  voltages.  High  operat¬ 
ing  voltage  prevents  full  oxidation  of  the  fuel,  which  re¬ 
sults  in  low  electrical  current  per  mole  of  incoming  fuel. 
This  region  is  characterized  by  high  thermodynamic  effi¬ 
ciency  but  low  Faradaic  efficiency.  At  peak  electrical  effi¬ 
ciency,  a  compromise  between  operating  voltage  and  fuel 
utilization  is  made.  Operation  at  peak  efficiency  results  in 
considerable  fuel  slip.  For  the  steam  reforming  case,  the 
exhaust  gas  at  peak  efficiency  contains  3.4vol%  H2,  while 
POX  exhausts  3.3  vol%  H2  (or  0.14  moles  H2  exhausted  per 
mole  CH4  feed),  and  syngas  results  in  4.0vol%  H2  in  the 
exhaust. 

The  thermodynamic  limit  of  fuel  utilization  for  any  op¬ 
erating  voltage  Vcell  is  reached  when  Vcell  =  Ev,  i.e.  current 
will  not  flow  spontaneously  unless  the  operating  voltage  is 
below  the  reversible  voltage.  Once  Vceii  =  Ev,  no  current  is 
drawn  and  gas  composition  in  the  anode  channel  does  not 
change.  Clearly,  one  cannot  arbitrarily  increase  the  cell  op¬ 
erating  voltage  to  increase  electrical  efficiency.  An  impor¬ 
tant  finding  is  that,  high  open  circuit  voltages,  while  de¬ 
sirable,  should  not  be  used  to  infer  high  electrical  effici¬ 
ency. 

4.  Conclusions 

In  the  spirit  of  Carnot  efficiency  for  heat  engines,  we  have 
developed  a  numerical  model  which  predicts  the  theoretical 


electrical  efficiency  for  SOFCs.  Since  Gibbs  minimization 
is  used  to  calculate  gas  composition,  any  fuel  mixture  may 
be  considered  without  a  priori  knowledge  of  the  chemical 
reactions  involved.  The  model  assumes  that  the  cell  length 
is  sufficiently  long  for  gas  composition  in  the  anode  channel 
to  reach  equilibrium  with  the  Nernst  potential  (Er  =  Vcell)  at 
the  outlet. 

When  a  second  assumption  holds,  i.e.  transport  of  ions 
through  the  electrolyte  is  the  rate  limiting  step,  the  model  is 
able  to  predict  conditions  along  the  length  of  a  SOFC  anode 
channel.  These  conditions  include: 

•  current  density  and  power  profiles; 

•  gas  phase  mole  fraction  profiles; 

•  fuel  slip  as  a  function  of  cell  length,  operating  voltage,  or 
changing  inlet  flow; 

•  minimum  fuel  cell  length  for  a  given  electrical  efficiency. 

SOFCs  promise  exceptional  efficiency  when  operated  appro¬ 
priately.  Nonetheless,  our  results  show  that  one  cannot  ar¬ 
bitrarily  increase  the  operating  voltage  to  increase  electrical 
efficiency.  Peak  electrical  efficiency  is  a  compromise  between 
fuel  utilization  and  operating  voltage,  which  are  antagonistic. 
We  also  find  that  high  open  circuit  voltage  is  not  necessarily 
a  good  indicator  of  high  electrical  efficiency.  We  conclude 
that  it  is  possible  to  exceed  60%  electrical  efficiency  in  in¬ 
ternally  reforming  SOFCs.  We  also  conclude  that,  all  things 
considered,  POX  may  provide  adequate  efficiency  for  many 
SOFC  applications. 
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